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We study the time reversal and space inversion symmetry properties of those transfer matrices 
mostly used in the calculation of energy spectra and transport-process quantities. We determine 
the unitary transformation relating transfer matrices. We consider the Kohn-Liittinger model for a 
quasi-2D system and show that even though the system studied in the (4 x 4) scheme satisfies all the 
symmetry requirements, the (2 x 2) subspaces do not fulfill such constrains, except in the F point 
of the Brillouin Zone. We find new exchange properties between the (2 x 2) subspace quantities. 
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I. INTRODUCTION 

In the last years a large number of theoretical studies of multilaycred semiconductor structures, using different 
schemes of k ■ p multiband Hamiltonians, have been published 1-17 . Besides the extensive application of these models, 
the symmetry problem in general, and the specific characteristics in particular, such as the time reversal invariance 
in the Kramer (2 x 2) subspaces of the Kohn-Liittinger (KL) model and the parity of heavy- (hh) and light-hole (Ih) 
states, have been only partially considered in a small number of references 1 ' 12-15 . The increasing interest in studying 
transport properties within the fe ■ p approximation puts forward the need of a simple procedure to obtain scattering 
amplitudes from Envelope Function Approximation (EFA) models, and to establish clear requirements to preserve 
the fundamental physical symmetries. Many physical properties, such as tunnelling and quantum coherence, depend 
strongly on the symmetries preserved by the system-particle interaction 18 . Recently, remarkable efforts have been 
done to study the breaking of discrete symmetries 12,13 ' 19-23 . However some issues remain to be studied. The lack of 
a theoretical analysis of whether the time reversal symmetry and the space inversion symmetry in the Kramer (2 x 2) 
subspaces are fulfilled or not, shows the need of a detailed and comprehensive analysis of symmetry problems in the 
KL and other EFA models. This is the aim of the present paper. 

An important quantity, where the original Hamiltonian symmetries can be checked out in a simple way is the 
transfer matrix (TM). This object is being used with increasing frequency to study transport properties of different 
type of systems and within quite distinct approaches. Basically, two types of transfer matrices are being used in the 
literature. One of these transfer matrices, called here of the first kind, relates a field F(z) and its derivative F'(z) at 
any two points z\ and Z2 of the scatterer system 26-28 . Another transfer matrix, called here of the second kind, connects 
state vectors i&(z) that are expressed in the propagating modes representation 27,29-34 . Because of its simple functional 
relation with the scattering amplitudes, the latter matrices were mainly used in the scattering approaches to deal 
with transport processes. The rigorous fulfillment of symmetry requirements, such as time reversal invariance (TRI), 
spin rotation (SR), space inversion invariance (SII) and flux conservation (FC) principle, were amply discussed in the 
literature in connection with these transfer matrices, both for 1-D one-channel and for 3-D multichannel systems 3 . 
Transfer matrices of the first kind or their associated transfer matrices (which relate F(z) and a certain linear form 
G{z) = v(z)F'(z) + u(z)F(z), at any two points Zi, and Z2), were mainly used in connection with the solution of 
EFA models 35,36 . The matrices u(z) and v(z) appear as coefficients in the matrix equation of motion 36,37 . To obtain 
the transfer matrix in the EFA one needs to build up the vector F(z), called the envelope (see Refs. 13,27,35,38-42 
and references therein). TM of the first kind have been widely used to study the hole spectrum in III-V, II- VI and 
IV semiconductors 24 ' 25,43-47 . 

Since the multiband KL model in the k ■ p approximation is amply used in the current solid state physics, and 
the understanding of the transport properties in heterostructure systems become of great interest, it is convenient to 
establish a clear bridge between the first and the second kind of transfer matrices. To this purpose the analysis of the 
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symmetry requirements on the transfer matrices of the first kind and their relation with the scattering amplitudes and 
transfer matrices of the second kind, is very much called for. We apply our results to study with clear advantages and 
transparency the multichannel transport of heavy and light holes through III-V semiconductor heterostructures 48,49 . 
However the results that we present here are to some extent also applicable to any EFA model. 

Assuming the existence of iV-coupled-differential-equations solutions, we define in Section II the transfer matrices 
of the first and second kind and determine the unitary transformation between them. In Section III, we derive the 
mathematical conditions that the transfer matrices of the first and second kind have to satisfy to fulfill the physical 
symmetries, and we also recall those for the TM of the second kind that are known. We focus our attention on the 
time reversal invariance, space inversion symmetry, and the flux conservation principle. In Section IV we consider the 
(4 x 4) KL Hamiltonian and derive an explicit representation of the unitary transformation determined in Section II. 
Next, we discuss the consequences of the various symmetry requirements obtained in Section II on the Hamiltonians 
and related quantities of the KL (2 x 2) subspaces. Finally we present some concluding remarks and comments. 



II. TRANSFER MATRICES: AN OUTLINE OF BASIC DEFINITIONS 



In this Section, we will outline the well known transfer matrix definitions. Based on these definitions, we will 
establish the relation between the transfer matrices of the first and second kind. 

Solving the system of equations associated to an EFA Hamiltonian, one obtains a set of linearly independent 
eigenvectors. It is often convenient to consider, as will be the case here, an orthogonal set of eigenvectors. With these 
eigenvectors we can build the envelope function 



2N 



where N is the dimension of the system of coupled differential equations, and the vector 



F(z) 
F\z) 



Alternatively, we can write the eigenvectors in the propagating modes representation as 

/*(*) = 



92j 



V 9Nj 

and define, in the plane-wave basis, the state vectors 

*(*) = 



a ip (z) 
by (z) 



(1) 



(2) 



(3) 



(4) 



where <p (z) and tp (z) represent the right and left moving modes (when qj is real). Both of them are A^-dimensional 

vectors built in terms of the right- and left-travelling wave functions <Pi (z) and ¥?j (z) in the i-th channel or 
propagating mode. 

It is easy to see that, writing <&(z) as 



*(*) 



ip 

o Xp 



(5) 



where a and b arc N dimensional vectors, and 



<P (z) 



( 4x (z) ^ 

V 2 (z) 







(0) 



\ o o 

a simple and obvious unitary transformation can be established between ^f(z) and <f?(z). Since these vectors are used 
to define the transfer matrices of the first and second kind, one can determine such unitary transformation, i.e. 

*(z)=Af®(z) (7) 

which is fundamental for the purpose of this paper. The specific structure of A/" depends on the particular Hamiltonian. 
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A. Transfer matrices of the first kind 

Using the function ^f(z) it is common to define the transfer matrix of the first kind M fd such that 

*{z 2 )=M fd (z 2 ,z 1 )*{z 1 ). (8) 
A slightly different definition for this type of transfer matrix is the associated transfer matrix T 36 defined by 

^(z 2 )=T(z 2 ,z 1 )^(z 1 ) (9) 

where 



- F(z) - 




G(z) 





F(z) 

v{z)F'{z) +u(z)F(z) 



(10) 



These two matrices are the first kind of matrices we use. 



B. Transfer matrices of the second kind 

Another TM that has been widely used in the literature that keeps a simple relation with the scattering amplitudes 
is the matrix defined by 31 



*(*a) = M sv (z 2 ,z 1 )^(z 1 ). 



(11) 



As for the first kind of transfer matrices, there is also a slightly different definition of the transfer matrix, usually 
called the coefficients transfer matrix defined by 



d =QM U 



(12) 



Based on the previous definitions and the transformation (7) it is clear that both types of transfer matrices are related 
by the system-dependent unitary transformation 



M sv (z R ,z h )=Af l M ld {z n ,z- L )M. 



(13) 



This is an important result. Based in the previous unitary transformation it is easy to obtain simple functional 
relations with the scattering amplitudes. The transformation matrix A/", say 



Af = 



9i 92 

di9i d-292 



92N 
d-2N92N 



(14) 



can be obtained when each linearly independent (LI) solution is written as a (N x 1) spinor, with no coordinate 
dependence (represented here by g^), times a plane wave. By dj we denote the coefficient of z in the exponent of the 
plane waves. In Section IV we will consider a particular KL Hamiltonian and we will obtain specific expressions for 
M fd and M sv and the unitary transformation in (13). 

III. GENERAL FC, TRI AND SII REQUIREMENTS ON M/ d AND M s „ 

Let us now recall and deduce the constrains imposed on the transfer matrices by the fulfillment of flux conservation 
and the physical symmetries of time reversal and space inversion. For the sake of brevity, we report the most important 
relations and quote just some of the necessary and well known results. 



A. Flux Conservation 



To obtain the flux conservation requirements we need to write the current density 

ih 



2m* 



(15) 
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in terms of any of the previously introduced 2iV-dimcnsional spinors. The method of derivation follows the usual 
calculation of the particle current in quantum mechanics. Wc thus have 



F'G — G~F 



(16) 



It is found useful to emphasize that in the KL model this formula reduces to j(z) = 21tl\[F'^ vF } — 2F'u^F. In 



other interesting cases as the Schrodingcr equation and the 1-D one-channel and for 3-D multichannel systems 
reduces to the known expression j(z) = F^F' — F'^F. 

When flux is conserved, which means j(z-i) = j( z i) at any two points, we immediately obtain the identity 



30-33 



R}(z)SyR(z) = M\ d (z,z )RHzo)-Z y R(zo)M fd {z,z ), 



(17) 



where 



R(z) 



and 



u(z) v(z) 



tin Om 



Henceforth O n /J n is the corresponding (N x N) null/identity matrix. In a completely similar way, one obtains 
the FC requirement for other transfer matrices. For the transfer matrices of the second kind the current conservation 
implies 



dc* — cd* = ba* — ab* 
which, using the transfer matrix definition, leads to the well known relation 

= Ml v (z,z )i: z M sv (z,zo), 



(18) 



(19) 



where S 2 = er 2 Jjy, is the enlarged (2N x 2N) Pauli a z matrix (see Appendix B). This condition has been explicitly 
deduced in various references 30-32 . 



B. Time Reversal Invariance 



TRI implies that the envelope function fulfills the relation 

KCF{z) = F{z) , 



(20) 



where the usual notation for the complex-conjugation operator C is used, and K is a model-dependent matrix 27,50 . 
Using the definition (8) in equation (20) we obtain 



M /d (z,z ) = SM} d (z,z )S- 1 



(21) 



with 



K O, 
O, K 

T, M sv and Q fulfill similar identities. In the case of the matrix T, we obtain 

T(z,z ) = £T*(z,z )5r 1 . 



(22) 



Concerning the transfer matrix of the second kind, it is well known that the TRI condition for systems with spin- 
independent interactions takes the form 



M sv (z, z Q ) = X x M* v (z, «b)£- \ 



(23) 



while for systems with spin-dependent interactions, this relation changes slightly. The matrix Y, x is given in the 
Appendix B. In the particular case of systems of spin 1/2 (sec Ref.32), TRI implies the condition 



M sv (z,z ) = KM* sv (z,Zo)K~ 



(24) 
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C. Space Inversion 



Let us now assume that the system possesses a space inversion symmetry with respect to the origin. It is convenient 
to emphasize that the space inversion we are interested here is the transformation that changes only the sign to the 
coordinate perpendicular to the interfaces of the Q2D system, leaving the signs of the in-plane wave vector components 
k x and k y unchanged. If ft denotes such space inversion operator, the space inversion invariancc implies 



SjF(z) = psF(- 



(25) 



where p = ±1 and s is a model-dependent matrix. 

As the derivatives change the parity of the functions, the action of the space inversion operator on the bi-field *f?(z) 
becomes 



with 



S = 



s 4 
4 -s 



Therefore, the SII requirement on the transfer matrix Mfd is 

M fd (z,z ) = S- 1 M fd (~z,-z )S. 
Similar identities can be obtained for all the transfer matrices. In particular, for M sv we obtain: 

M sv (z,z ) = SJvM sv (-z,-z )S sv . 
In the Appendix B, the explicit expression for S sv is given. 



(26) 
(27) 

(28) 
(29) 



IV. TM PROPERTIES AND SYMMETRY REQUIREMENTS IN THE KOHN-LUTTINGER 

HAMILTONIAN 

The first part of this section is devoted to obtain an explicit form for A/" in the case of the (4 x 4) KL model. In 
the second part we will focus our attention on the symmetries of the (2 x 2) subspaces of this model. For the sake of 
completeness we shall briefly recall some important and well-known results within the (4 x 4) KL model 1,12-15,51 ' 52 
and we will describe explicitly the consequences of its TRI and SII symmetries on the (2 x 2) subspaces. These 
consequences have not received sufficient attention, however, because of their current interest concerning the clastic 
and assisted transmission studied within the (2 x 2) subspaces, deserve further clarification. We shall discuss these 
implications for the cases in which the wave vector k (defined in the interface planes) is either equal or different from 
0. 



A. Derivation of the transformation matrix Af in the (4 x 4) KL model 



Dealing with the (4 x 4) KL model Hamiltonian it is usual to block-diagonalize it 5 and to work then with (2 x 2) 
instead of (4 x 4) Hamiltonians. In this way, the mathematical difficulties are highly simplified. This method remains, 
nevertheless, a very useful tool to understand many of the intriguing physical properties of hh and Ih valence bands 
near the band edge, and we use it for the derivation of the transformations (7) and (13) in the (4 x 4) space as our 
main objective. 

In the block-diagonalizing procedure, a unitary transformation U performs the splitting of the original (4 x 4) 
Hamiltonian H into two (2 x 2) blocks, which arc labelled "up"(u) and "low"(Z) 

UHU ] UF(z) = 

where the blocks are given by 

A lK 2 + B 2 k 2 + V(z) C xy -iD xy k z 

C xy + iD xy k z A 2 k 2 +B 1 k 2 z +V(z) 



7~tu 2 
2 Hi 



E 



(30) 



7~tu = 



6 



H 



i = 



A 



2K 



B x kl + V{z) C xy -iD xy k z 



C xy + iD xy k z A lK 2 + B 2 k 2 + V(z) 



(32) 



For the definitions of A\. A2, B\, B 2 , C xy , D xy and k 2 see Appendix A. Thus, for each block and for each slab of 
material (with certain set of phenomcnological parameters 53 ), we have an eigenvalue problem 



{h u ,i-EI 2 }f u ,i{z) = 2 . 



Solving this equation we have 



(33) 



F{z) = U ] 




Fi(z) 
F 2 (z) 



(34) 



The spinor components are given in the Appendix A. Here aj(Pj) stand for the linear combination coefficients 
of the LI solutions times the corresponding normalization constant in the "up" ("low") subspaces, respectively. We 
remark that only two of the values of the momenta qj arc LI. They arc obtained from the zeros of the secular 
fourth-order polynomial determinant of (33), which is the same for both u and I Hamiltonians and then, the energy 
eigenvalues of the subspace hole states are the same with different spatial trends. This is perhaps the most important 
consequence of the block-diagonalization (B-D) process from (30). 

When considering the multichanncl-multiband transport of hh and Ih through, say, a III-V semiconductor 
heterostructures 48,49 , one should keep in mind the order of the basis components to assign correctly the transmission 
amplitudes. Assuming 

F{r) = F{z) ■ u {r) Q^(^+ y S v ) ^ 

we shall choose for the periodic part of the Bloch function Uq{t) = (m, u 2 , u 3 , 1*4) = (| |, |) , 1 1, — i) , 1 1, |) , 1 1, — |)), 
in agreement with the order hh +3 / 2l lh-\/ 2 , lh + \/ 2 , hh_ 3 / 2 given in Ref.5. Accordingly, the wave function F{z) must be 
supplemented by identifying correctly each value of the momenta qj to an hh or Ih state. To perform this identification 
we refer to the levels of an infinite quantum well (iQW) at n = limit, given by previous calculations 6,11,13 . Considering 
a parabolic approximation to the dispersion law for the hh and Ih bands in this limit, the energy can be cast as 



E h h = (71 - 272) q 3 
Em = (71 + 272)9?, 



(35) 



which unambiguously relates the uncoupled levels hh(lh) to (73(91), respectively. 

Before going forward, notice that both transfer matrices Mfd and T relate states in the quasi-particle representation 
characterized by a certain set of quantum numbers and describe the evolution along z of these modes with no mention 
to their propagation direction. Nevertheless, as suggested above we can alternatively describe the system using 



the propagation modes representation 



:!2 



the representation where right and left moving states are resolved into 



separated components. This idea in combination with the mentioned adjustments lead us to express the wave function 
(34) and its derivative in terms of an eight-component state vector (5) 



tf (z) = W 



G (G)* 
(*G)* 



if (z) O4 
O4 £ (z) 



(36) 



where 



W 



v o 4 
o 4 V 



V 






e-"> 
000 e-'* 



being 



- iarctan(— -) — arctan( 
2 I k y 



1l?,k x k y 



I2{k 2 



and 



= -|arctan(-) + arctan( 72(fcg _ fcg) ) 



7 



The parameters r\ and <f> are taken in the form required by the unitary transformation 5 U for block-diagonalization 
(30). Hereafter 71, 72 and 73 are the Liittinger parameters of the layer. We have defined 



G = 



513 .911 .941 343 

323 .921 .931 333 

"323 —321 .931 333 

"313 -.911 .941 343 



and O = 



q 3 

qi 

31 

q 3 



In (36) the (4x1) vectors a(b) contain the corresponding a,j(bj) linear combination coefficients of the LI solutions 
times the corresponding normalization constant of the configuration space in the "up" ("low") subspaces, respectively. 
Defining 



Af = W 



we have 



*(*) 



•A/" 11 



W*ATl2« 



G (G)» 
iGQ (iG)* 



(p (z) O4 
O4 V (*) 



(37) 



JV*(z). 



(38) 



This relation is essential to achieve the purpose posted at the beginning of this section. Thus, from (38), the expression 
(8) that relates ^f(z) at any two points Z\ and Z2 of the scatterer system can be written as 



{M $(z 2 )} 2 = M fd (z 2 , z 2 ) {Af &(z 1 )} 1 
Hence the state vectors (5) in these regions are connected as 

#(* a ) =Ar^ 1 M fd (z 2 ,z 1 )Ar i ^(z 1 ), 

therefore 

M sv (z 2 , Zl ) = Afj 1 Mf d (z 2 , Zl ) Mi . 



(39) 
(40) 
(41) 



The convenience of this (8 x 8) matrix identity is likely found in the clear advantages of using both types of transfer 
matrices for the description of tunneling and related properties of particles moving through multilayered structures. 
This result becomes a profitable platform 48,49 to circumvent the usual assumptions and restrictions that the current 
descriptions 11,56-60 face in describing the hole tunnelling in hctcrostructurcs. 

We emphasize several points here: Attention must be paid to orto-normalization of the spinor components 3^, 
content in the (8 x 8) transformation matrix Af, to prevent prompt failures during numerical simulation of quantum 
transport phenomena of heavy and light holes through heterostructures. For systems well-described within the EFA 
the orto-normalization procedure is not a trivial question and differs from routine process, basically, by the presence 
of linear k z terms in (31,32). A detailed study on that subject is published elsewhere 54,55 . 



B. TRI and SII symmetries for the k/0 case in the (2 x 2) KL model 

Let us analyze the implications of the TRI and SII symmetries (in the (4 x 4) KL model) on some physical quantities 
defined in the (2 x 2) subspaces. We shall start with the TRI symmetry. To exhibit this implications on the time 
reversal operator operator we will use its representation reported by Pasquarello et al. in Ref.l. Before block- 
diagonalizing the time reversal and space inversion operators, we need to transform them from the Pasquarello's basis 
representation to the Broido and Sham basis representation 5 . This transformation is done by means of 



-10 

10 

i 

-i 



(42) 



In the following, a x ,a y ,a z are the three Pauli matrices and the coordinate origin will be placed in the inversion 
symmetry point [z = 0). 
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1. Time Reversal Invariance 



The (4 x 4) time reversal operator in the Broido and Sham basis is given by 



K^aC — RK p R C R 



R 



2 tr y 

(Ty 2 



R T C 



2 <r x 

-(T x 02 



c. 



(43) 



Hereafter subindcxes BS and PB stand for Broido-Sham and Pasquarello-Bassani, respectively. After applying the 
transformation 5 



U = 



V2 



o e - *" -e" 1 o 
o e~ ir > e ir > o 
Q-i<i, o o 



(44) 



the time reversal operator T BS becomes 



T RS = URK m R T U T C 



02 cr x 

-cr x 2 



C. 



(45) 



As a consequence the TRI requirement implies for the (2 x 2) subspace Hamiltonians the fulfillment of the following 
condition 



cr x H u (z)a x = Hi(z) 
o- x Hj (z)cr x = H u (z) 



(46) 



It is then clear that the TRI of the KL Hamiltonian in the (4 x 4) space does not imply the TRI of the Hamiltonians 
in the (2 x 2) subspaccs. It is easy to show also that 



and hence 



here: 



<r*Ki(z)=Fi,u(z). 

r x M* ,(«)r x - Mi, u (z) 



(47) 
(48) 



(Tx 02 
02 <Tx 



As can be seen from previous relations, the considered physical quantities of the subspaces u and I should be related 
to each other in order to satisfy the TRI requirement in the (4 x 4) space. 



2. Space Inversion 



It is convenient to emphasize that the space inversion symmetry we are dealing with is the one which changes 
only the sign to the coordinate perpendicular to the interfaces of the Q2D system. The usual space inversion 1 makes 
{x,y,z} — > {— x, — y,— z} then, when taking 2D-Fouricr transform in the [x,y] plane, change the signs of k x and k y 
also. 

As in the case of time reversal, we need to know first the space inversion operator (25) for the (4 x 4) KL space in 
the Broido and Sham basis 5 . This is obtained from 



*5bs — RS PB R 



R 



er z 02 
02 er 2 



R 



I2 O2 
O2 I2 



(49) 



It is easy to verify that the (4 x 1) wave function F(z) defined in (34), under this operator action, satisfies the 
symmetry condition (25). It can be seen also that the time reversal and space inversion operators satisfy {T BS , S BS } = 0, 
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i.e. they anti-commutc. This leads to the existence of an additional symmetry operator S BS T BS that changes a wave 
function for a given n with a definite parity into the same function of k but with opposite parity 1 ' 13 . The invariance 
of the wave function F(z) in (34) under space inversion operation (25) implies 







-Fi(-z)' 






F 2 (z)_ 


= P 


F 2 (-z)_ 




F 2 (z)_ 



Since 



2 cr x 
&x 2 



the space inversion invariance of the KL Hamiltonian leads us to 



2 cr x 
<?x 2 



H u (-z) 2 
2 Hi(-z) 



2 cr x 

&x 02 



H u (z) 

o 2 



Therefore 



o 2 

ni(z) 



H u (z) = cr x Hi(-z)a- x 

Hi(z) = <T X H U {-Z)(T X 

For the transfer matrices of the first kind and for the envelope wave functions we have 

M uJ (z) = T x M Lu (-z)T x , 



(50) 



(51) 



(52) 



(53) 



<r x F u> i{-z) = pF lt u (z). 



(54) 



Here again the transfer matrices and wave functions, appearing in the first and second members of equations (53) 
and (54), belong to different subspaces. For k ^ the space inversion symmetry in the (4 x 4) space does not imply 
space inversion symmetry in the (2 x 2) spaces. Hence, the wave functions in both subspaces does not have a definite 
parity. 

The mixing of hh with Ih states increases as we go away from the center of the Brillouin Zone, i.e., as k differs from 
0. The square modulus of the wave function (T> u , i) was evaluated by using the representation for the TM of first kind 
for the hole state. Taking advantage of (54) the following useful relation is obtained: 



Vi(E,k,z) =V u (E,k,-z). 



(55) 



where T> u i are defined as 



t r 



T> u> i(E,k,z) = M"£ (z, z<^T u , i(z ) M u l2 {z,z )T' u l {z ) 



(56) 



where M^'^, z D ) arc blocks of the TM of first kind M"' l , respectively. 



C. TRI and SII symmetries for the k = case in the (2 x 2) KL model 

At k = the hole states decouple in the case of no external field or no strain, given independent series of levels 
(35) for the infinite quantum well boundary conditions. Both lri u and Hi are diagonal and thus invariant under time 
reversal and space inversion operations. The time reversal operator in this case is <t z C, which commute with the 
Hamiltonian, and the wave functions satisfy 

cr z F Utl (-z)=pF U! i(z), (57) 

then the hh and Ih components of the vectors F u _i(z) posses definite but opposite parities. At the T point the up 
and low subspaces satisfy FC, TRI and SII and belong to the symplectic universality class 32 , as the original (4 x 4) 
space of solutions. 
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v. 



CONCLUSIONS 



We establish the relation between the TM of the first kind (frequently used to determine the energy spectra) 
and the TM of the second kind (widely used in the scattering approach) in general and for the specific case of KL 
Hamiltonian model. This relation allows us to take the advantage of using the appropriate transfer matrix according 
to the circumstances one is dealing with. It is possible to avoid the usual assumptions and restrictions that the current 
descriptions face in describing tunnelling of holes 11 ' 56-60 . 

We deduce new symmetry requirements on the TM of the first kind M fd and SII requirements on TM of the second 
kind M sv for N > 1 case of actual interest in quantum-transport problem. Any EFA model having a master equation 
like that of the References 36,37 are expected to have the same or similar properties. In fact, a great part of the 
analysis carried out in sections II and III can be extended to those problems with minor changes. 

We made clear that the symmetry analysis in the reduced (2 x 2) subspaces of the KL model demand attention to 
prevent erroneous statements 11 . We show that for k ^ the time reversal and space inversion invariance requirements 
of the (4 x 4) KL Hamiltonian does not imply the invariance of the 'up" and "low" Hamiltonians and related objects, 
separately. However, to fulfill the whole Hamiltonian invariance they have to transform one to the other according to 
the corresponding expressions obtained here. It has been also deduced that for k = 0, the "up" and "low" subspaces 
recover continuously the time reversal and space inversion symmetries, as expected. 

Besides the basic interest of this topic, the relations obtained here arc useful and convenient in controlling bother- 
some numerical calculations, as well as in studying the relevant quantities in the multichanncl-multiband tunnelling 
of holes within a more realistic approach. Expressions of section III were verified numerically within computer uncer- 
tainty. 
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APPENDIX A: PARAMETERS OF THE (2 x 2) H u AND Hi KL HAMILTONIANS 




(Al) 



The non-ortonormalized spinor components gij 



.9ii = 

.913 = 
.922 = 
.924 = 
.931 = 
333 = 
342 = 
.944 = 



A 2 K 2 + B x q\ + V(z) y 
A 2 k 2 + B iq 2 + V(z) 
-(C xy + iD xy q 2 ) 
~{C xy + iD xy q 4 ) 
A\K 2 + B 2 qj + V(z) 
A\n 2 + B 2 q 2 + V(z) 
-{C xy + iD xy q 2 ) 
-(C xy + iD xy q 4 ) J 



(A2) 
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from 7t u and Hi peculiarities follow 



#9 = 


-91 1 

"2 1 


t^ 


— <7s 

TO 


.7 11 


jli 


(7"| Q = 


fil /I 

14 


.922 = 
.924 = 


.921 . 

* / 

523 


.931 = 


532 


.933 = 


534 


542 = 


041 


.944 = 


543 , 



(A3) 



being 

gij,g 3 j real 1 
52j ,54j complex J ' 

The Liittingcr parameters 71,72,73 typify each layer of the hcterostructure. 



APPENDIX B: TRANSFORMATION SYMMETRY MATRICES 

Given the transformation between the matrices Mfd(z, Zo) and M sv (z, zq), it is easy to see that the FC constrain 
on M sv (z, zq) (Sec. Ill) implies that 

S z = (Af^ 1 JfdAf. (Bl) 

However, it is worth noticing that, whenever the basis of linearly independent KL Hamiltonian eigenvectors, is not 
an orthogonal basis, we will have 

it* (2) _ u ( z ) —v(z) 
fd ~ [ v(z) O4 _ ' 

For TRI and SII requirements on the TM M sv (z, z ) (Sec. Ill) we used the matrices T, x and S sv , respectively, which 
read 

S x =Af- 1 Y,- 1 Af* 7 Ssv = ATSAf- 1 . (B2) 
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